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ABSTRACT 

Due to collective processes small solid particles 
exhibit very strong resonances in absorption and 
scattering in the spectral regions in which the 

real part of the dielectric constant, is negative. 

The precise positions depend on the particle shape, 
size and on the nature and amount of coatings it 
might have. A comprehensive discussion of these 
resonances is given in the Rayleigh approximation. 

It is shown that the observed bands at 2200 A 
and in the 10 y region are most probably due to such 
collective excitations. Some other processes are 
also considered and are shown to be relatively inef- 
ficient and inadequate. The following specific con- 
clusions are drawn: (a) the 2200 A interstellar band 

is very likely due to graphite particles; (b) these 
graphite particles should be very small, approximately 
spherical, and should have no coating whatsoever; (c) 
the identification of circumstellar and interstellar 
silicates from the observations in the 10 y region 
does not seem to be correct; (d) very valuable infor- 
mation about the shape of the circumstellar and inter- 
stellar dust particles can be obtained directly from 
observations; (e) narrow band polarization measurements 
in the spectral regions of these bands will be very 
helpful in determining the shape of the particles. 


I. INTRODUCTION 

Recent OAO-2 observations of the interstellar extinction 
curve have been discussed by Bless and Savage (1972; see also 
this volume) . The observations show large intrinsic varia- 
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tions . Bless and Savage have concluded that a multi-component 
model of interstellar grains, such as I had proposed earlier 
(Gilra 1970, 1971), is needed to explain the observations. 

The most significant feature in the entire extinction curve 
is the ultraviolet bump at about A -1 = 4.6 y -1 . The observa- 
tions indicate that, in general, this feature is very strong, 
has a quite symmetrical profile, and is not "flat-topped". 

Also, the variations in the position of the extinction maximum 
are very small. In this paper our main concern is to identify 
the physical process and the substance which can give rise to 
such a feature at this wavelength (about 2170 A) . 

It seems, as will be shown below, that the most plausible 
process which can satisfactorily explain this band is that of 
plasma oscillations. This process was introduced in astronomy 
by Unsold (1964) to explain the diffuse interstellar band at 
4430 A. He also mentioned that the position of the resonance 
due to this process (plasma resonance) shifts to longer wave- 
lengths as the size of the particles is increased, but he did 
not explain this effect. It was pointed out by van de Hulst 
(1964) that the strength and position of plasma resonances de- 
pend on the shape of the particles. Apart from the work of 
Unsold (1964) and that of van de Hulst (1964) there is no dis- 
cussion of this process in the astronomical literature. Simi- 
lar resonances occurring in the infrared for some solids will 
be discussed. Observations also indicate a band in the 10 y 
region (Woolf and Ney 1969, Hackwell et ad_. 19 70) . 

It will be shown that the theoretical resonances occur due 
to "collective" processes. A comprehensive discussion of 
these resonances will be given in the small-size limit. In 
this limit they can be studied with the help of solutions of 
relatively simple problems in electrostatics. The following 
cases will be considered: homogeneous ellipsoids, coated el- 

lipsoids, and ellipsoids with anisotropic optical properties. 

A very brief discussion of the effect of particle-size on the 
position of the resonances will be given. The physical pro- 
cesses which give rise to these resonances will be mentioned. 

Next, the astronomical observations in the ultraviolet and 
infrared will be compared with the theoretical results. Some 
other physical processes will also be considered and the sug- 
gested identifications will be critically examined. Finally, 
some observationally verifiable predictions will be made. 

II. THEORY 

a) Small Homogeneous Ellipsoidal Particles 

Let us discuss absorption and scattering by homogeneous el- 
lipsoidal particles in the Rayleigh approximation. For a field 
applied along one of the main axes (j = 1,2,3) of the ellipsoid. 
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the polarizability a- is a tensor and is written as (van de 
Hulst 1957, 1964) J 


cu = V 

4tt Lj (£-e m ) + e m 


(la) 


where V is the volume of the particle, e its dielectric con- 
stant, e m is the dielectric constant of the medium, the Lj de- 
pend on the ratios of the axes, and the values 4 ttl^ are known 
as the depolarization factors. The Lj satisfy thecondition 
0 1 Lj <. 1# with the property j Lj = 1. If m = n - ik is the 
complex refractive index, we have e = m 2 = ej - ie 2 = 
n - k 2 - i2nk. If the particle is in free space, e m = 1 , 
and equation (la) reduces to 


a . . V £-1 

3 W Lj(e-l) + 1 


(lb) 


Equations (la) and (lb) apply to the case of long elliptical 
cylinders also (van de Hulst 1957, p. 71) . 

The scattering cross-section, and the absorption cross- 

section, C^, are given by 


and 



( 2 ) 

(3) 


where Re means "the real part of". Putting the value of ou 
from (la) in (2) and (3), we get J 

r j = 8 tt 3 v 2 ( £ l -e m) 2 + e 2 2 

sea 3 A 1 *" {Lj£^+e m (1-Lj) } 2 + (L j£ ; 2 ) 2 ' 


and 

c 3 _ 2 ti_ „ £ m £ 2 . 

abs A {LjE^ejjjd-Lj ) P + (Lj e 2 ) 2 ' 5 

Here we have assumed that is real; we shall also assume 
that its value does not change with wavelength. 

From equations (4) and (5) we see that there is a resonance 
in the absorption and scattering cross-sections at the wave- 
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lengths where 


Lj£^ + (1— Lj) — 0 , (6) 

and Lj £2 < 1* This is because, as pointed out by van de Hulst 
(1964;, the polarizability of the particles is very large as 
can be seen from equation (la) . As we see from equation ( 6 ) 
there is one resonance corresponding to each value of Lj . We 
therefore have three resonances for an ellipsoid, two for a 
spheroid and one for a sphere. These resonances occur where 
£]_, the real part of the dielectric constant, is negative and 
are confined roughly within the maximum of £ 2 / corresponding 
to Lj = 0, and the maximum of (e^ 2 + £2 2 ) -1 / corresponding to 
Lj =1. This is the case of a thin disc. The location of the 
particle in a medium of dielectric constant £ m or in free 
space ( £jn = 1 ) does not affect the position of these bounda- 
ries. However, for intermediate values of Lj , the resonances 
in the case of free space are at 

£ L = - { (1/L j ) - 1} , (7a) 

and in the case of a medium at 

e l = " e m { C 1/L j ) - 1} . (7b) 

When discussing real solids we shall see that these conditions 
imply a shift in the position of the resonances to lower ener- 
gies (or longer wavelengths) as the particles are placed in a 
medium of dielectric constant £ m because of the way £^ varies 
with wavelength. As an example we take the case of a spheri- 
cal particle. L 1 = L 2 = L 3 = 1/3, and in free space the reso- 
nance is at £^ = - 2 , whereas if the particle is in a medium of 
dielectric constant £ m , the resonance occurs at £^ - 2 £ m , 

that is, at a higher negative value of £ 3 . 

We may remark in passing that C sca has a minimum which oc- 
curs on the higher energy side of the resonance. The condi- 
tions from equation (4) are £^ = £ m and £2 = 0. These are the 
same conditions which give a minimum in the reflectivity, and 
the frequency at which these are satisfied in the infrared is 
known as the Christiansen frequency. 

b) Size Dependence 

As the size of spherical particles is increased the resonance 
shifts to longer wavelengths (Unsold 1964) . Similar results 
are obtained for spherical graphite particles (see, e.g. , 

Gilra 1971) . For very small sized graphite particles the ex- 
tinction maximum is at about 2100 A and as the mean size is 
increased the feature becomes broader and the maximum shifts 
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to longer wavelengths. We have to use the Mie theory to ex- 
plain this (see, e.g., Doyle and Agarwal 1965, Ruppin and 
Englman 1968). A detailed discussion of size dependence will 
be given elsewhere (Gilra 1972) . We just note that as the 
first departures from the Rayleigh approximation become appar- 
ent and we consider only the electric dipole term, the reso- 
nance occurs at a higher negative value of £]_ given by 

£1 = -2 - 11 x 2 , (8) 

5 

where x = 2i\a/\ and a is the radius. This result can be rea- 
dily obtained from an equation given by van de Hulst (1964, 
p. 28), though he does not seem to have noticed this. 

As the mean size becomes quite large the resonance becomes 
broad and shallow as can be seen in the extinction curves of 
SiC (Gilra 1971, Figure 2). The broad shallow feature short- 
ward of about 1800 A seen in that figure is due to such a 
resonance. Once again this resonance is to the high energy 
side of the maximum of ^2* 

c) Effect of Coatings 

Let us consider an ellipsoid of dielectric constant £i hav- 
ing a coating in the shape of a confocal ellipsoid and a die- 
lectric constant £ 0 . Let this compound ellipsoid be placed in 
free space. Following Stepin (1965) and Bilboul (1969) we can 
derive an expression for the polarizability, a, for such an 
ellipsoid. For a field applied along one of the main axes the 
polarizability is given by 

v i 

a = V_ (e 0 -l){L i (e i -e 0 ) + £ 0 > + -y-( ej-e 0 ) iL ( l-e Q ) +£ Q ) } 

47T {(e i -e 0 )L i +e 0 Hl+L(e 0 -l) } + L(l-L) (e r e 0 ) (e Q -l)' 

V 

where Vj_ is the volume of the inner ellipsoid, V is the total 
volume, and 4uLj_ and 4 ttL are the depolarization factors of the 
inner ellipsoid and the total ellipsoid respectively. For the 
sake of clarity we have suppressed the direction dependence. 

If Vj_/V = 0 , or Vi/V =1, or £j_ = e Q , or £ 0 = 1/ this equation 
reduces to equation (lb) for the case of homogeneous ellip- 
soids. If = 1, we get an expression for the polarizability 
of a hollow shell whose inner and outer surfaces are confocal 
ellipsoids . 

In the degenerate cases, thin discs, spheres and thin 
needles, we have Lj_ = L. For coated spheres L^ = L = 1/3 and 
the .result is identical to that obtained for the small size 
case from the rigorous solution (Guttler, 1952, van de Hulst 
1957, p. 74). For a thin disc, or a highly flattened oblate 
spheroid, L^ = L and there are two values of L, L = 0 and 
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L = 1. The results are quite interesting. 
For L = 0 , 


a = 


V 

471 


{ (1-Xi) (e -1) + 6, -1) } 


V, 


V 


V 


or. 


a={— (e 0 -l) + — (£■; -1) } , 

4tt ° 4tt 


and for L = 1, 


V r „ V ±> e Q 1 ^i e i”l- 


a = - {(1 —) — 
4tt V £ 


V £ ■ 


or. 


a = { 


Vq £q ^ Vi £ i — 1 


+ 


1 '-x 


4tt e 0 4 tt e- 


L } , 


( 10 ) 


( 11 ) 


where V is the volume of the shell and Vj_ is the volume of the 
inner ellipsoid. On comparing equations (10) and (11) with 
equation (lb) we note that the polarizability of such a com- 
pound disc (or, for limiting values of L) is equal to the sum 
of the polarizability of two separate thin discs having the 
volumes V Q and , and with dielectric constants £ 0 and 
respectively . 

Scattering and absorption cross-sections in the Rayleigh ap- 
proximation for compound ellipsoids can be calculated by using 
the value of polarizability from equation (9) in equations (2) 
and (3) respectively. The case of a compound thin disc is 
simple; we can use equations (10) and (11) in equations (2) 
and (3) . For such a thin disc as in the case of polarizabili- 
ty, the absorption cross-section also is the sum of the absorp- 
tion cross-sections of two separate thin discs having the 
volumes V Q and Vj_ , and dielectric constants £ Q and £i respec- 
tively. For other cases the final expression for C a b s and 
C SC a are somewhat complicated if both £j_ and £ 0 are complex 
but they become relatively simple if £ 0 is real. Before we 
write the expressions, let us compare the results for the case 
of a thin needle, or a highly elongated spheroid (L = L-^ ; 

L = 0 and L = 1/2) , with those for a coated infinite circular 
cylinder in the small radius limit. 

From the rigorous solution of Kerker and Matijevic (1961) 
for perpendicular incidence, I have obtained the following 
dominant terms (up to x 2 ) in the approximation that the total 
radius is much smaller than the wavelength both outside and 
inside the cylinder. For the TE mode (E 1 symmetry axis) the 
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dominant term is 


— 


_ 1TTX 


(Sq- 1) (Sj+So) + q 2 (sj-Sp) (1 +Sq) 

„ (Ci+tQ) (l+e 0 ) + q 2 (£i-e 0 ) (e 0 -l) . 


= nrx 


say , 


(12a) 


and for the TM mode (E || symmetry axis) it is 


b = H2L 

o A 


11TX 


-' { (1-q 2 ) (e 0 -l) + q 2 (ei-1) } 


B, say, 


(12b) 


where x = 2irr/X, r is the total radius, q = rj[/r, and r^ is 
the inner radius, and as before, e Q and are dielectric con- 
stants of the shell and core respectively. These reduce in 
the appropriate limits to the case of a thin homogeneous infin- 
ite cylinder. C e xt/ the extinction cross-section per unit 
length and C sca f the scattering cross-section per unit length 
for the two cases are given by 


TE mode : 


TM mode: 


^ext = 
^sca = 
^ext - 
C sca " 


2r 

2r 

2r 

2r 


2_ 

x 

2 

x 

2 

x 

2 

x 


Re (2a x ) 

* 2 I a i I 

Re (b Q ) 



r 


(13) 


(14) 


Let us now compare the cross-sections for a thin infinite 
circular cylinder with those of a small highly elongated pro- 
late spheroid. The TE mode corresponds to L = 1/2 and the TM 
mode to L = 0 for the highly elongated prolate spheroid. In 
terms of A and B, as defined in equations (12a) and (12b) , the 
polarizabilities for a small highly elongated prolate spheroid 
are , 


for L = 1/2, a = (V/4 tt) 2A, (15) 

and for L = 0, a = (V/4 tt)B, (16) 

where we have made use of the fact that in this case V^/V = q 2 . 

If £ is the major axis and r the semi -minor axis of the whole 
spheroid, V = 2/3 iTr 2 £. The absorption and scattering cross- 
sections can now be calculated with the help of equations (3) 
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and (2) respectively. We note that for absorbing particles in 
this approximation C ex j- = C a fr s and we can compare the two in a 
straightforward way. 

So that a direct comparison is possible the cross-sections 
for the two cases, thin infinite circular cylinder and small 
highly elongated prolate spheroid, are presented in a tabular 
form. We note that for the former case the tabulated cross- 
sections are per unit length . 


Thin Infinite 
Circular Cylinder 

Small Highly Elongated 
Prolate Spheroid 

TE mode: 

C ext = 2£*-irr 2 *2 • [Re(iA) ] 

A 

Csca = 1T 2 r * (27rr/A)3 * t l A l 2] 

L = 1/2: 

Cgbs = ^-*jffr 2 Jl«2 • [Re(iA) ] 
C S ca = ^y 71 ^ 2 • ( 27Tr/X) 4 • [ | A | 2 ] 

TM mode: 

L = 0: 

c ext = ^*7rr 2 • [Re (iB) ] 

9tt 9 

C abs " X*3 7rr2£ * [Re(iB)] 

C sca “ ^•(27,r/X) ! -[|B| ! l 

C sca “ ^H 2 -( 2ur/X>‘.[|B| ! ] 


^ext anc ^ ^abs f° r the two cases are the same; the only dif- 
ference is that in the corresponding volumes. For the scat- 
tering cross-sections the terms in the square brackets, and 
these are the more important terms, are the same because of 
the depolarization factor, but the other terms are different. 
For the infinite cylinder case we have a X -3 dependence, 
whereas for the highly elongated prolate spheroid the wave- 
length dependence is X -If . (Foi all ellipsoids in general in 
this approximation the scattering cross-section has a X -lt de- 
pendence.) The same is true for the corresponding homogeneous 
cases (cf., van de Hulst 1957, p. 316; Kerker 1969, p. 266). 

We now write expressions for C a bs and C sca for the compound 
ellipsoid on using equation (9) in equations (3) and (2) if e Q 
is real. « 




( 18 ) 
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where we have substituted (e-^-i£ 2 ) for e i/ and 


R = L ± (e 0 -l) + ^ {e o - L(e 0 -1)> , 

S = (1-Li) C£ 0 " 1 > " " L ( £ o _1)} = ^o" 1 ) - R ' 

T = + L(e 0 -1)} + (e 0 .-l)-jL L(l-L) , 

and U = (1-Li> {.1+L (e Q -l) } - (e 0 -l) Yi. L (1-L) = {l+L(e Q -l)} - T. 


We shall now make the further assumption that £ 0 is inde- 
pendent of wavelength. Let us investigate the nature of the 
resonances.. For the limiting cases L = L-^ and, L = 0 and 
L= 1 and the results can be easily obtained from equations 
CIO) and (11) .. For C a ^ p the coatings do not have any effect 
on the position or the strength of the resonances but for C sca 
the shell is also contributing and there is a slight change. 
For intermediate values of L, 0 < L < 1, the resonances occur 
at a value of £]_ given by 


e l 



L + L(e 0 -1) ' 

Li + L(e 0 -l){Li + ^L(l-L)} 


(19) 


It follows from equation (19) that if there is no coating the 
resonance occurs at = -(-1 + 1/Li) as in the case of a 
homogeneous ellipsoid. It can be shown that as the relative 
amount of coating is increased (keeping V fixed and decreasing 
Vi) the resonance occurs at higher negative values of eq. Fi- 
nally, if Vi/V -*• 0, the resonance occurs at ej = -e Q (-l+l/L^) , 
corresponding to the case of a homogeneous ellipsoid in a 
medium of dielectric constant £ a . Thus whatever the amount of 
coating the resonance is always between = -(-1 + 1/Li) and 
£1 = — e 0 ( — 1 + 1/Li) • This result is similar to that for sur- 
face plasmons in the case of coated thin foils (Raether 1967) 
and is more general. (These results apply to the case of a 
coated thin infinite cylinder also.) As noted earlier, the 
boundaries of the resonances do not change because of coatings 
We also note that the minimum in C sca does not shift from 
ei = 1 to £]_ = £ m because the shell is also contributing to 
the scattering; the shift is actually to the lower energy side 
The result we shall be making use of in this paper is that 
even a slight coating makes a substantial change in the posi- 
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tion of the resonance. 

Although the formula for polarizability, equation (9), is 
implicit in the work of Stepin (1965) and Bilboul (1969) , most 
of the discussion in this section is believed to be new. 

d) Small Ellipsoids with Anisotropic Optical Properties 

Since graphite is a suggested interstellar particle and 
since its optical properties are highly anisotropic (that is, 
e is a tensor) , let us discuss the nature of these resonances 
for such particles in the small size limit. The polarizabili- 
ties can be calculated from the work of Jones (1945; see also 
van de Hulst 1957, p. 73). We consider an oblate spheroid 
whose axes coincide with the principal directions of the di- 
electric constant tensor. (The results can be extended to the 
more general cases.) If e ± and e„ are the dielectric con- 
stants when fi, the electric vector of the incident radiation, 
is perpendicular and parallel to the symmetry axis respective- 
ly, and 4 ttLj_ and 4 irL t | are the respective depolarization fac- 
tors, the corresponding polarizabilities, a j_ and a ( |, are given 
by (Wickramasinghe 1967, Greenberg 1968) 


V e x -1 

aj - 4 tt L_ L (e 1 -l) + 1 


a„ 


_V_ £||-1 

4 ir L |( (£,, -1) + 1 


( 20 ) 


where V is the volume of the particle. These equations are 
similar to equation (lb). As before, we can calculate the 
scattering and absorption cross-sections from equations (2) 
and (3) for the two cases separately. 

Let us write e ± = e ± , i “ i £ x ,2 and £ n = e i.i ,1 ” i £ n , 2 • As 
in the case of a homogeneous oblate spheroid, there are two 
resonances, one for fi perpendicular to the symmetry axis (Elc) 
and the other for fi parallel to the symmetry axis (fillc) , 
assuming, of course, that both £j_ x and e | ( ^ 1 become negative 
in some energy ranges. There are* several interesting aspects. 
Even for a spherical particle, which is one extreme case of an 
oblate spheroid, there are two resonances, one for file at 
e x r 1 = -2, and the other for fine at e,, x = -2. As the parti- 
cle becomes oblate the resonance for ijj_c shifts to higher 
negative values of x and for fillc to lower negative values 
of e n . Finally, for the limiting case of a thin disc 
(platelet or highly flattened oblate spheroid) , the resonances 
are at the maximum of e lr2 and the maximum of (ejf /1 + e,,^) -1 
respectively. 
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e) Physical Processes 

So far we have discussed the nature of these resonances in 
terms of £, the dielectric constant. We noted that in spite 
of the possible changes in the shape, amount of coating, size 
and optical anisotropy, the resonances a^e confined within the 
maximum of £2 the maximum of (e£ + £2) _1 * If all these 

parameters are specified the resonance occurs at a fixed nega- 
tive value of , and hence for a given solid, at a fixed 
energy or wavelength. In other words, the excitations are 
quantized. The question now arises: what are the energy 

ranges in which is negative for the solids of astrophysical 
interest? 

Depending on the type of solid, becomes negative, in 
general, in two energy ranges (Kittel 1966) : one in the infra 

red due to lattice vibrations and the other in the ultraviolet 
The behavior can be quite complicated and in each of these 
energy ranges there can be more than one such region. In the 
infrared the frequency corresponding to the maximum of £2 is 
called toip , the transverse optical phonon frequency, and the 
frequency at which £^->0 and £ 2+0 is called w L , the longitudi- 
nal optical phonon frequency. Wip is less than u)l. In the 
ultraviolet the frequency at which £p->0 and £ 2^0 is called oop, 
the plasma frequency. There is no frequency corresponding to 
(jjp, but there is a maximum in £2 due to the transitions of 
bound electrons on the lower energy side of cop. 

We now very briefly discuss the physical processes involved 
in the excitation of the resonances we have been discussing in 
this paper. If infrared photons are incident on a diatomic 
crystal of ions ±e, the positive and negative ions are dis- 
placed with opposite signs. Quantized vibrations that follow 
are designated as "material polaritons" (Englman and Ruppin 
1968a, b; Ruppin and Englman 1968) . In the ultraviolet the 
electrons take the place of negative ions . The electron gas 
starts moving as a whole with respect to the positive ion 
background. This is known as plasma oscillation and the quan- 
tum of this excitation is termed a plasmon (Kittel 1966) . 

Both these processes can be characterized as collective excita 
tions . Thus the resonances discussed earlier in this paper 
can be identified as occurring due to collective processes. 

Two qualifications should be made, however. The peak in z2 in 
the ultraviolet is, as noted earlier, due to the transitions 
of bound electrons. Therefore the resonances occurring at 
this energy for some shapes are not due to a collective pro- 
cess. In the case of gold and silver in the visual and near 
ultraviolet, for example, there is what is known as the free 
carrier plasma frequency, but there is no corresponding peak 
in £2 on the lower energy side of it. The discussion of the 
resonances has to be modified somewhat to take this fact into 
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consideration for such cases. 

The discussion of the physical processes involved has been, 
of necessity, very simplified. For details reference should 
be made to papers cited in this section (see also, Kawabata 
and Kubo 1966, Steinmann 1968). 

Let us now summarize all the theoretical discussion on 
these resonances. They are due to collective processes, with 
the exceptions noted above, and are confined within io T and 
in the infrared and within the maximum of £2 an< ^ W P i- n the 
ultraviolet. The precise position depends on the shape, size 
and amount of coating. We can invert this argument. From an 
observed feature, if it can be established that these process- 
es are involved, we can determine, at least in principle, the 
shape and the size of the particles and also how much coating 
these particles have. 1 

III. COMPARISON WITH OBSERVATIONS 
a) The 2200 A Band 

As earlier calculations indicated there is an extinction 
band for small-sized spherical graphite particles at about 
2200 A. This has been attributed to transition of tt electrons 
to the conduction band in the graphite crystal (Wickramasinghe 
1967, Stecher 1969) . This identification is incorrect because 
this band occurs where e ± , 1 is - 2 , and is therefore due to 
plasma oscillations in spherical graphite particles. The di- 
electric constant of graphite in the basal plane, £j_ , is given 
by Taft and Phillip (1965, and private communication) and is 
plotted in Figure 1 for energies between about 2.5 ev and 9.0 
ev. The plasma frequency, Wp, is at about 1800 A (7.0 ev) . 

The maximum in £3 is at about 2800 A. For spherical particles 
the resonance is at £]_ = -2 which is at 2205 A. For an oblate 
spheroidal particle with the ratio of minor to major axis = 0.8 
the resonance shifts to 2255 A. Similarly, if we coat a 
spherical graphite particle of 50 A radius with a shell of ice 
about 5 A thick, the resonance shifts to 2255 A. 

. In Figure 2 is plotted the absorption cross-sections using 
equations ( 20 ) for graphite particles for two extreme shapes, 
a sphere and a thin disc, normalized such that they are unity 
at the maximum. If C^g and C^k s are the absorption cross- 
sections on using e x and e n respectively, the quantity plotted 
is ( 2 C^k s + C^bgJ/S, as is the case if the particles are ran- 
domly oriented. Values of e i( for graphite were taken from 
Tosatti and Bassani (1970) . There are no strong resonances in 

c abs* 

The highest negative value of e'^ is -0.85 which occurs at 
12.7 ev. In this approximation C^s is greater than C^j-^ be- 
tween about 1300 and 850 A, the exact wavelengths depending on 
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Figure 7 . — The complex dielectric constant ol graphite In the 
b(U>al plane [*ee text lor dlAcvutlon) . 


the shape. The net effect is to increase somewhat the total 
C a bs i- n this wavelength range from what is obtained from C^b s . 
But at these wavelengths the Rayleigh approximation is not 
valid and we cannot say what the contribution is for larger 
sizes because the required theory is not available. It may be 
remarked that in C^bs there is a much stronger resonance in 
the far ultraviolet than the one at 2070 A. This resonance 
occurs at about 800 A and is also due to the process of plasma 
oscillations . 

The peak in the small size limit as calculated with equa- 
tions (20) occurs at 2080 A for spheres; however, e x f x is 
equal to -2 at 2205 A. The shift to shorter wavelengths could 
occur for two reasons. From the formulas we see that the reso- 
nance occurs where Z2 also satisfies the condition, Z2 < 1* 

It occurs at shorter wavelengths if this condition is not 
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satisfied (see, Kawabata and Kubo 1966). The second factor is 
that Cgkg is also contributing and it is increasing towards 
shorter wavelengths, the net effect being to push the reso- 
nance to shorter wavelengths. The second factor is not impor- 
tant in this case, however. 

Two other computational details should be noted about such 
resonances in general. If one does not anticipate where the 
resonances will be and does not compute cross-sections for a 
sufficient number of wavelengths, one is likely to get mis- 
leading results. Secondly, one should use the dielectric con- 
stant and not the refractive index to get physical and compu- 
tational insight about these resonances. 

The peak for spherical particles occurs at 2080 A and as 
the oblateness increases it shifts to longer wavelengths. 
Finally, for thin discs it occurs at 2780 A. Also plotted in 
Figure 2 is the observed interstellar extinction curve in the 
region 2860-1840 A for cr Sco (Bless and Savage 1972) using two 
different base lines for its normalization. The broader one 
(large dashes) is obtained if the base line is parallel to the 
abscissa and starts at the furthest point observed in the 
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infrared. The narrower one (small dashes) is obtained if the 
base line is a smooth curve joining 2.6 y -1 and 6.5 y -1 in the 
extinction curve of o Sco. The width of the observed feature 
depends critically on the choice of the base line. The peak 
in the observed curve is at 2170 A and the feature is quite 
narrow. 

It is obvious from Figure 2 that the observations put very 
stringent limits on the shape of the particles. We can defi- 
nitely rule out graphite platelets (thin discs) as being re- 
sponsible for the observed feature. 

My calculations for spherical particles, using the Mie the- 
ory with £j_ , show that for a radius of 0.005 y the peak is at 
2070 A, for 0.010 y at 2100 A, and for 0.020 y at 2200 A. If 
we now compare the theoretical curves taking all the factors 
into consideration with the observed curves, we can conclude 
that the suggested interstellar graphite particles should be 
fairly spherical, quite small (mean radius ~ 0.015 y) and have 
almost no coatings whatsoever. The minimum in the observed 
curves between about 1800-1500 A is also consistent with the 
graphite hypothesis (Gilra 1971) . 

If we take the values of e from Scouler (1969) for Mg 2 Si, 
the calculations show a strong feature at about 2200 A for 
spherical Mg 2 Si particles due to plasma oscillations. The 
list could be made substantially bigger but it is not clear if 
the particles of other type will be formed anywhere. 

As we have noted, the position and strength of the reso- 
nance depend critically on the precise values of e (A) . How- 
ever, in the discussion so far we have tacitly assumed that we 
can use the bulk dielectric constants for particles in inter- 
stellar space. This may not be exactly true: for sizes ~ 100 

A, £ depends on size due to quantum effects which tend to 
broaden the plasma resonance (Kawabata and Kubo 1966) ; we do 
not know what effect the temperature (the temperature of inter- 
stellar particles is ~20°K) has on £; the interstellar parti- 
cles may not be identical to bulk crystals in their crystal- 
line structure. Surface "irregularities" also affect the 
nature of plasma resonances (Steinmann 1968) . We still do not 
know how to calculate the extinction cross-section for parti- 
cles with anisotropic optical properties if the size is arbi- 
trary. The effects of these and perhaps some other related 
problems should be investigated. 

Another suggestion to explain the 2200 A band is that it is 
caused by an absorption edge. Huffman and Staap (1971) ex- 
perimentally obtained the optical constants for one silicate, 
(MgFe)SiC> 3 , up to 12.4 ev. They suggested that the observed 
feature could be explained with a gaussian size distribution 
having a mean size 0.06 y. In Figure 3 are plotted extinction 
(E) , scattering (S) and absorption (A) cross-sections for 
spherical particles of this silicate for two mean radii. 
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0.06 y and 0.07 y. The suffixes refer to the mean radii in 
microns. The curve E # q 6 is the same which according to Huff- 
man and Staap explains the observed feature. 

At an absorption edge the imaginary part of the refractive 
index (or the imaginary part of the dielectric constant) in- 
creases rapidly. At this energy the absorption cross-section 
increases very rapidly and the scattering cross-section 
decreases very rapidly, as can be seen from the curves . For 
very small-sized particles the extinction is entirely due to 
absorption. The relative contribution from scattering in- 
creases as the mean size is increased. For a narrow size- 
range a "band" appears in the scattering cross-section due to 
the rapid decline on the high energy side (as already noted) 
and a rapid rise on the low energy side. With only a 15% in- 
crease in the mean size, the rapid rise on the low energy side 
occurs at much lower energies making the scattering "band" ex- 
tremely broad. The resulting "band" in the extinction cross- 
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section is "flat-topped", asymmetric and extremely broad and 
bears no resemblance in shape, strength or position with the 
observed feature. Therefore it seems unlikely that such a 
well-defined observed feature can be caused by such a physical 
"process". Wickramasinghe and Nandy (1971b) have also noted 
some of these points independently. 

Another physical process could be considered. Some solids 
either due to an intrinsic property or due to an impurity or 
due to "physical damage" exhibit absorption bands in their 
bulk optical properties. For example, neutron irradiated 
quartz shows a band at about 2200 A (Billington and Crawford, 
Jr. 1961) . Let us investigate if such particles can explain 
the observed feature. (Dr. G. H. Herbig also mentioned this 
possibility during the discussion following my paper.) Van de 
Hulst (1957, p. 191) has already discussed such a theoretical 
problem. Let us refer to Fiqure 3. The feature in the absorp- 
tion cross-section at 3.8 y _r is due to such a process. Only 
in the very small-size limit does it show up in the extinction 
cross-section. As the size increases the minimum in the scat- 
tering cross section makes it invisible. In fact, for the 
mean size 0.07 y there is a minimum in the extinction cross- 
section at this wavelength, as already discussed by van de 
Hulst. It could still be argued, however, that only very 
small-sized quartz particles may be responsible for the ob- 
served feature. However, there is a stronger band in neutron 
irradiated quartz at 1650 A which is not seen in the observed 
extinction curves. Also, since quartz has an absorption edge 
at somewhat higher energies (the scanner #1 on board the OAO 
has a quartz window and that is why no observations are possi- 
ble with it shortward of about 1800 A) , the observed extinc- 
tion curve should show a very rapid increase at higher ener- 
gies. The observations of o Sco and p Oph down to 1100 A 
(Bless and Savage 1972) which show strong 2200 A features be- 
lie this expectation. Thus it does not seem likely that such 
particles are responsible for the observed band. 

Let us summarize the discussion on the identification of 


* One remark should be made about my model of interstellar 
grains (Gilra 1971) . In this model silicates are the chief 
contributor to the far ultraviolet rise. This conclusion was 
based on the assumption, which seemed most expedient at the 
time, that the optical constants are wavelength independent. 
As discussed here the optical constants are now available for 
one silicate. The calculations show that the suggested exis- 
tence of silicates to explain the far ultraviolet rise may be 
somewhat uncertain now (see Bless and Savage (1972) for other 
possibilities) . 
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the 2200 A band. Three different physical processes were con- 
sidered: plasma oscillations, absorption edges, and absorp- 

tion bands in bulk crystals. It has been shown that the last 
two are relatively inefficient and inadequate. In these cases 
the behavior of scattering is almost opposite to that of ab- 
sorption in the region of interest (Figure 3) . As the size is 
increased beyond a certain critical value the extinction band 
becomes exceedingly broad and highly asymmetric in one case, 
and invisible in the other. On the other hand, the band due 
to plasma oscillations arises because there is almost identi- 
cal structure in the bands in both scattering and absorption. 
Whatever changes we may make in the particle shape (Figure 2) , 
size (Figures 1 and 2 in Gilra 1971) and coatings (as dis- 
cussed earlier) the extinction band is always well defined and 
is very strong in the small-size limit. It can be concluded, 
therefore, that of all the suggestions made so far the process 
of plasma oscillations is the most plausible process to give 
rise to a feature of the type observed. It is correct, how- 
ever, that quite stringent limits were placed on the geometri- 
cal and surface properties of the suggested interstellar gra- 
phite particles because the observed feature is quite narrow. 
Some theoretical problems were mentioned earlier and they 
should be investigated. It may be noted that the minimum in 
the observed curves between about 1800-1500 A is also consis- 
tent with the graphite hypothesis (Gilra 1971) . 

b) The Bands in Infrared 

It has been suggested (Friedemann 1969a, b; Gilman 1969) 
that particles of silicon carbide can form in the atmospheres 
of stars in which the carbon to oxygen ratio is >1. In the 
spectra of many of these stars the violet region is strongly 
depressed, the quasi-discontinuity setting in at about 4400 A. 
It has been suggested (Stephenson and Ross 1970) that this 
phenomenon is "among the most outstanding of all opacity ef- 
fects to be seen in stars of any spectral type." Silicon car- 
bide has the fundamental absorption edge between 2.2 and 3.0 
ev, the precise position depending on the polytype, impurities 
and temperature. Gilra and Code (1971, 1972) have suggested 
that circumstellar silicon carbide grains are responsible for 
the observed violet depression. Their theoretical results are 
in very good accord with the observations. The model of 
interstellar grains I have proposed (Gilra 1971) has silicon 
carbide as a major constituent. We are therefore interested 
in investigating the nature of resonances in small silicon 
carbide particles in the infrared so that the observations 
could confirm or disprove the suggested existence of such cir- 
cumstellar and interstellar particles. 

The fundamental absorption band in SiC is at 12.6 y 
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(Spitzer <at al . 1959a, b) , that is, w T , the transverse optical 
phonon frequency, corresponds to this wavelength. w L , the 
longitudinal optical phonon frequency, corresponds to 10.3 y. 
(The optical properties of cubic and hexagonal SiC are quite 
similar in this region and we shall consider only the cubic 
form.) My Mie calculations for spherical particles had shown 
(Gilra 1971) a weak band at but this was completely dwarfed 
by a very strong absorption feature at "'10.75 y where eg = -2. 
I had suggested that the feature will become broad if the par- 
ticles are not exactly spherical. At that time I was not 
aware that experimental results existed (Pultz and Hertl 1966) 
which confirm this prediction. The theory has been discussed 
in several papers by Englman and Ruppin (Englman and Ruppin 
196 8a, b; Ruppin and Englman 1968) . 

Figure 4 shows the absorption cross-section for cubic SiC 



Figune 4. — Shape, dependence o& ne*o nance* in SiC in the 10-13 y 
region Uee text ion di-6 ciu^ion) . 704 in the uppen night 

connen should nead 79 4. 
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particles of three shapes, thin disc, sphere and needles, 
using equation (5) with e m = 1 . Random orientation has been 
assumed. For a thin disc the resonances are at 970 cm -1 
(= 10.3 y = w L ) and at 794 cm -1 (= 12.6 y = w T ) . As the 
oblateness decreases the resonances move towards each other 
and finally for the spherical case there is only one resonance 
at 932 cm -1 (= 10.73 y; = -2) . As the particle becomes 
prolate, the resonances move away from each other and finally 
for the case of a needle they are at 794 cm -1 (= 12.6 y = w T ) 
and 948 cm -1 (= 10.55 y; = -1). If the particles are 
"highly irregular" there will be a broad feature between about 
10.3 y and 12.6 y such as has been observed experimentally by 
Pultz and Hertl (1966) . The results for needles are also in 
very good agreement with their experimental results for fibres 
We have no a priori way of determining the shape of the 
suggested circumstellar and interstellar silicon carbide par- 
ticles. There should be a strong absorption band (or bands) 
between 10 y and 13 y in the interstellar extinction curves of 
the highly reddened stars . This feature should appear in emis 
sion in the stars which show strong violet opacity. However, 
in the latter case this predicted emission will be somewhat 
"filled in" by strong absorption in this wavelength region due 
to SiC 2 molecules (which give rise to the Merrill-Sanford 
bands in the late N-type carbon stars) and high spectroscopic 
resolution will be needed to separate the two effects (Gilra 
and Code 1972) . 


c) The 10 y Band and Silicates 

A broad emission feature centered at about 10 y has been 
seen in the spectra of oxygen-rich late type variables (Woolf 
and Ney 1969) . Perhaps a similar feature at 10 y has been 
seen in absorption in the observations of the galactic center 
(Hackwell et_ ad. 19 70) . These .features have been attributed 
by the authors to circumstellar and interstellar silicates re- 
spectively. Since the dielectric constant for silicates is 
not yet available in this spectral region, let us investigate 
the theoretical basis for these identifications. 

If n, the real part of the refractive index, is very close 
to one, and k, the imaginary part, is very much smaller than 
one, then equation (5) reduces to (van de Hulst 1^57) with 
£ m ~ 1 

Cabs = ( 2 tt/X) V* 2k 

and since the absorption coefficient A = (4irk)/A, we get 

Cabs = 


V* A . 


(21) 
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Equation (21) gives the absorption cross-section in the Ray- 
leigh-Gans approximation. This is the equation used by Woolf 
and Ney (1969) and Hackwell et. al . (1970) , following the work 

of Gaustad (196 3) , in identifying the silicates. Donn et_ al . 
(1970) have agreed with this assumption. But this is not cor- 
rect because the conditions n - 1 and k << 1 are not satisfied 
for silicates in this region. The structure in ei and £2 
should be somewhat similar to that for quartz (experimental 
absorptance results (Lyon 1963) show this very well) and that 
means very high values of £2 and £ 1 (both positive and nega- 
tive for £]_) . The dielectric constant for quartz is given by 
Spitzer and Kleinman (1961) and theoretical calculations 
(Knacke 1968, for spherical particles, Martin 1971, for sphe- 
rical and cylindrical particles, and my own for ellipsoidal 
particles) using equation (5) show no match at all with the 
values given by Gaustad (1963) using equation (21). (Most of 
the values in these spectral regions given by Gaustad (1963) 
are incorrect for other substances as well, for the same rea- 
sons.) There are strong resonances (as in the case of SiC) 
which depend on the shape of the particles. The same should 
be the case for silicates. For silicates Wickramasinghe and 
Nandy (19 71a) have taken the value of k from Gaustad (196 3) 
and assumed a wavelength independent value of n (= 1.66) 
throughout this spectral region. This assumption is also not 
correct. We thus conclude that the identification of circum- 
stellar and interstellar silicates from observations in the 10 
y region is based on a wrong theoretical premise and inade- 
quate experimental data and therefore seems to be premature. 

There is one way, however, in which the experimental ab- 
sorption coefficient can be compared directly with the obser- 
vations, and that is if we assume that the distribution of 
shapes and sizes of the particles and their spatial arrange- 
ment are identical in the two cases. It is not clear if this 
is the case. The best way to resolve the controversy is to 
get experimental values of £ and to obtain observations at 
high spectroscopic resolution to see the structure within the 
broad feature. Then the theoretical curves using the values 
of £ obtained experimentally can be directly compared with the 
observations . 

Some of the points in this section have been independently 
discussed by Martin (1971) . 

d) Polarization 

From equation (5) we see that absorption, and therefore 
thermal emission, depend on the shape and orientation of par- 
ticles. If the particles have any anisotropy, either in shape 
or in optical properties, and if they are suitably aligned, 
then the 2200 A band and the bands in the infrared should be 
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polarized. There should be a polarity reversal at some wave- 
length in these spectral regions . 

If the observed interstellar polarization in the visual is 
due to the same agent which is responsible for the 2200 A 
band, then this band must be very strongly polarized, as 
Stecher and Donn (1965) have already pointed out for graphite. 
Recent balloon-borne wide band polarimetry of the reddened 
star C Oph by the University of Arizona astronomers (Coyne et_ 
al . 1971, and remarks at the 135th meeting of the American 
Astronomical Society) indicates that the polarization contin- 
ues to decrease in the region 3000-2000 A. If these observa- 
tions are correct and if the observations of £ Oph can be con- 
sidered as representative of other stars as well, then we can 
conclude that the suggested interstellar graphite grains are 
not aligned and the polarization in the visual is not dime to 
graphite grains. 

If for the late-type variables the observed polarization in 
the visual is due to grains and the emission in the infrared 
is also due to the same grains and these grains are suitably 
aligned, then the emission feature(s) should be polarized with 
a polarity reversal at a fixed wavelength. However,, very nar- 
row band polarimetric observations will be required so that 
the polarizations with different polarity do not cancel each 
other. Dr. M. Dyck at the Kitt Peak National Observatory has 
started such an observational program. 

In the model of interstellar grains I have proposed (Gilra 
1971) SiC particles are responsible for the observed inter- 
stellar polarization in the visual. As discussed earlier SiC 
particles have resonances in the 10-13 y regions.. Therefore 
the suggested interstellar band(s) in the 10-13 y. region, should 
be polarized with a polarity reversal in between. 

It is obvious that narrow band polarization measurements in 
these spectral regions will give us important information about 
the shape, size and chemical composition of circumstellar and 
interstellar grains. 


IV. FINAL REMARKS 

Very strong resonances in absorption and scattering occur 
in the spectral regions in which ej_, the real part of the 
dielectric constant, is negative. The precise positions depend 
on the particle shape and size and on the nature and amount of 
coating. The nature of these resonances was discussed with 
the help of solutions of relatively simple problems in electro- 
statics and the discussion given in this paper is essentially 
complete in the small-size limit. In this limit they are al- 
ways confined within the maximum of e 2 on the low energy side 
and the maximum of (e.^ 2 + e2 Z )~ 1 on the high energy side. As 
the size becomes comparable to the wavelength, new resonances 



COLLECTIVE EXCITATIONS AND DUST PARTICLES IN SPACE 317 

appear on the low energy side of the peak in £2 (see Ruppin 
and Englman 196 8, Englman and Ruppin 196 8b).. These will be 
discussed elsewhere (Gilra 1972) . 

There are two distinct energy ranges in whirjh these reso- 
nances occur. One is in the infrared in the region of strong 
lattice bands between and w L . The other is in the ultra- 
violet between the maximum of £ 2 / which is due to the transi- 
tions of bound electrons, and wp. Broadly speaking, the pro- 
cess which gives rise to these resonances is that of collec- 
tive excitations . 

Taking all the factors into consideration and comparing the 
properties of these resonances with those of "bands" due to 
other processes, it can be concluded that these resonances are 
the strongest and, relatively speaking, the most permanent 
features in the whole domain of light scattering by small par- 
ticles . 

The discussion given in this paper can be applied to some 
problems in solid state physics, colloidal chemistry and 
astronomy. With these results it has been possible to explain 
(Gilra 1972) some experimental results in solid state physics 
and colloidal chemist ry which have not yet been accounted for. 

The conclusions concerning some astronomical problems are 
very interesting. It has been shown that the observed bands 
at 2200 A and in the 10 y region are most probably due to col- 
lective excitations. Graphite particles continue to be the 
best working hypothesis for the 2200 A band, even though strict 
limits have been placed on their geometrical and surface pro- 
perties because of the narrowness of the observed band. A very 
important result is that for the first time, at least in prin- 
ciple, it may be possible to determine the shape of the circum- 
stellar and interstellar particles directly from observations . 
Narrow band polarization measurements in the regions of these 
bands will be very helpful in achieving this goal. 

I am grateful to Dr. A. D. Code for his helpful advice and 
constant encouragement throughout the course of this work. 
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